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Abstract 

The problem of determining the acoustic field in an inviscid, isentropic 
fluid generated by a solid body whose surface executes prescribed vibrations 
is formulated and solved as a multiple scales perturbation problem, using the 
Mach number M based on the maximum surface velocity as the perturbation 
parameter. Following the idea of multiple scales, new "slow" spacial scales 
are introduced, which are defined as the usual physical spacial scale 
multiplied by powers of M. The governing nonlinear differential equations 
lead to a sequence of linear problems for the perturbation coefficient 
functions. However, it is shown that the higher order perturbation functions 
obtained in this manner will dominate the lower order solutions unless their 
dependence on the slow spacial scales is chosen in a certain manner. In 
particular, it is shown that the perturbation functions must satisfy an 
equation similar to Burgers’ equation, with a slow spacial scale playing the 
role of the time-like variable. The method is illustrated by a simple 
one-dimensional example, as well as by three different cases of a vibrating 
sphere. The results are compared with solutions obtained by purely numerical 
methods and some insights provided by the perturbation approach are discussed. 
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1. Introduction 


We wish to describe the acoustic field in an inviscid, isentropic fluid 
generated by a solid body whose surface executes prescribed vibrations. When 
the Mach number M, based on the maximum surface velocity is small, it is 
natural to consider a perturbation approach to the problem. Crow (1970) 
carefully analyzed a general class of problems of this type, including 
Lighthi 1 1 ’ s (1952, 1954) acoustical analogy approach, involving quadrupole 
sources, and Ribner’s (1962) approach, based on monopole sources. He showed 
how this class of problems could be treated as a singular perturbation problem 
by the method of inner and outer asymptotic expansions. Whitham (1956, 1974) 
examined this problem in the high frequency limit by using a perturbation 
approach similar to Lighthill’s (1949, 1961) method of strained coordinates. 

He essentially "strained" the time coordinate and obtained expressions for the 
acoustic field valid in the far field for high frequency excitations. The 
approach we shall use is not restricted by the high frequency assumption, 
although our results do reduce to Whitham’s results (to leading order in M) in 
the far field and at high frequencies. In addition, the analysis we shall 
present, being based on multiple spacial scales, provides insight into the 
effects of the nonlinearity on the form of the solution, as well as insights 
into the interaction of the value of the Mach number and the angular 
dependence of the solution (see section 9). 

In section 2 we formulate the mathematical problem we wish to solve and 
then outline the multiple scales approach in section 3. We apply our method 
to a simple one-dimensional example in section 4. In section 5 we formulate 
the problem for a vibrating sphere and then determine the leading terms in the 
multiple scales expansion in sections 6 and 7. In section 8 we apply our 
results to three different vibrating spheres and then discuss our results in 
section 9. 

2. Problem Formulation 

We wish to describe the acoustic field in an inviscid, isentropic fluid 
generated by a solid body whose surface executes prescribed vibrations. To 
begin, we let the origin of a Cartesian coordinate system (x^,X 2 ,Xg) be fixed 

at some convenient location inside the body. Then, in terms of these 


1 


coordinates, the equations of conservation of mass and momentum are 


( 2 . 1 ) dp/d t + ^«(pu) = 0 , 

(2.2) dp/dx j + p{d^iJdt + (u»^)u^j = 0, i=l,2,3 . 


In (2.1) and (2.2) p, u^ , p, and t represent the fluid density, the fluid 
velocity component in the positive x^ direction, the fluid pressure, and the 

-A. -X -i -pfc 

time coordinate, respectively. Also, u = u i. + u~i 9 + u„i„ and V = 


i (d/dx. ) + i 0 {d/dx 0 ) + i Ad/dx„) , where i. is a unit vector in the positive 

1 1 & U O J 


Xj direction. We also let p = k (P , where k and 7 are known constants, and 
define 


(2.3) c 2 (p) = dp/d p = k 7 (p)' 1 1 . 

We let pQ, U, and L be typical (constant) values for the density, 

velocity, and length scales, respectively, associated with the flow, and 
define nondimensional variables (denoted by a " A " above the quantity) by 


(2.4) 


A A A 2 

P = P/P 0 , “ = u/U, p = (P-P 0 )/(P 0 U ), 

Xj = Xj/L ( j =1 ,2,3) , t = t/(L/c Q ), 


where p Q = kp^ and c Q = c(p Q ). We now rewrite equations ( 2 . 1 ) and (2.2) in 

terms of these variables and then omit the " A " above the various quantities to 
obtain the relations 


(2.5) 


dp/dt + M^*(pu) = 0, 
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( 2 . 6 ) 


— 0 , i — 1,2,3 , 


dp/dx . + p 2 7 ^Mdu./dt + M 2 (u.V)u.j 

where M = U/c Q is the Mach number of the flow. 

We now seek approximate solutions to equations (2.5) and (2.6), subject 
to appropriate boundary and initial conditions, which will be formally valid 
for small values of M. 

3 A multiple scales perturbation solution of the basic equations 

Following the method of multiple scales (see e.g. Nayf eh( 1973) , Chapter 

^(k) 

5), we introduce the spacial scales (variables) x , k = 0, 1, 2, ... , 
related to x by 

(3.1) x (k) = M k x, k = 0, 1,2 


Thus x (0) = x, x (1) = Mx, x^ = M 2 x, etc. We then assume that u and p are 


functions of these new scales, as well as functions of t and M, i.e. p 

. ,M) . 


p(t,x (0) ,x (1) ,x (2) , . . . ,M) and u = u(t,x^°\x^ ,x^ 


(k) 


The method of multiple scales now treats all of the variables {x. } as 


independent variables and seeks to determine p and u as functions of these 

variables. In particular, for "small" values of M, we look for solutions in 
the form 


oo 


p = l + Mp (1) + M V 2) + . . . = , with P (0 = 1, 

(3.2) j =0 


oo 


u = u (0) + M^ (1) + M 2 U (2) + 


. = ^u ( j¥ . 
j=0 


3 



is independent of M, but, 


Here each of the coefficient functions p ^ and u^ 


in general, will depend upon t and the spacial scales xj k \ i.e. 

p (j> (t,i ( 0 ) ,i (1) ,...) and . 

To determine these coefficient functions, we substitute the expansions 
(3.2) into equations (2.5)-(2.6) and use the relation 


(3.3) V = V (0) + MV (1) + M 2 V (2) + ... , 

where V (k) s ^ (d/<?x| k) ) + i^d/dx^) + i^d/dx^) . 

We then collect coefficients of like powers of M on the left side of each 
equation and, hence, express the left side of each equation as a power series 
in M. We then equate to zero the coefficient of each power of M, since the 
right side of each equation is zero. In this way, we obtain the following 

system of equations satisfied by the coefficient functions p^^ and u^ ^ : 


(3.4a) dp {l) Id t + V (0) .u (0) = 0, 

(3.4b) dp^/dxj 0 ^ du|°)/dt = 0, i = 1, 2, 3; 

(3.5a) dp (k+1) /dt + V (0) .u (k) = F (k) , 

(3.5b) d/> (k+1) /dxj 0) + du| k) /dt = G[ k) , i = 1, 2, 3, k£l. 


(Equations (3.4) follow from the terms in equations (2.5)-(2.6) which are 
O(M) , while equations (3.5) follow from the terms in these equations which 

are 0(M k+1 ).) Here the functions and depend only upon with 

j<k+l and u^ with j<k. In particular, we find 

(3.6a) F (1) = -(V (0) .(p (1) u (0) ) +^ (1) .u (0) ), 
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(3.6b) 


g< 1) = -(<u (0) .v <0, )u‘ 0) + (2-'i)/ 1) auj 0) /at ♦ a/ 1) /<fc' 1) ). 


To solve equations (3.4), we set 


(3.7) 


1] = -dip/dt and ^ (0) = V ( °V 


Then equations (3.4b) are satisfied for any choice of <p, while equation 
(3.4a) yields the requirement that <f> must satisfy the linear wave equation 


(3.8) 


d 2 tp/d t 2 - A ( °V = 0, 


where (d/dxj 0 ^) 2 + id/dx^) 2 + (d/dx^) 2 is the usual Laplacian 

operator in the variables {x|^}. 

From the structure of equations (3.5), we see that these equations can 
(in principle) be solved recursively, starting with k=l . In particular, 

equations (3.5) are a system of linear equations for the unknowns p ^ ^ and 

u^^ . Consequently, we can express the solution to these equations as the 
superposition of a particular solution and a homogeneous solution to these 
equations. The homogeneous solution has the same form as the solution to 
equations (3.4) (see equations (3.7) and (3.8)). The particular solution 
will, of course, depend on the index k, as well as on the solutions for the 
lower order perturbation coefficients. In general, these particular solutions 
have the property (as we shall demonstrate explicitly in the following 

( 1 ) 

sections) that they tend to decay more slowly in magnitude as | x |->oo than p 


and u 


(0) 


Consequently, if this decay were to be left unchecked, the 


perturbation expansions (3.2) would become invalid as |x|^oo. However, as we 

(k+1) -*(k) 

shall show, it is possible to satisfy the requirement that p and u 

should decay at least as fast as p^^ and u^^ as |x|-»oo by properly choosing 

the dependence of p^+1) and on th e "slower" spacial scales x^ ^ , x , 

etc. Thus, this method will yield a perturbation expansion (3.2) which will 
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be uniformly valid as |x|->oo and, consequently, will correctly represent the 
far field behavior of the acoustical radiation. 

In the following section we shall demonstrate some of these ideas with a 
simple, one-dimensional example, and then proceed to a class of three- 
dimensional problems in sections 5-8. 

4. A Simple One-Dimensional Example 

As a simple example to illustrate the general ideas of our approach, we 
consider first a flow which varies in only one spacial coordinate, such as 
one-dimensional flow in a semi-infinite tube. We let this one spacial 
coordinate be denoted by x^ = x, so that equations (2.5)-(2.6) become 

(4.1) dp/d t + M(d/dx)(pu) = 0 , 

(4.2) dp/dx + p 2 ^M<?u/dt + M 2 u(du/dx)j = 0 . 


We shall assume that u is a specified function, say f(t), at x=0 for all t^O, 
and that ambient conditions prevail at t=0 for all 0<x<®, i.e. 

, . u| =f(t)for t^O, and u| = 0, p| = 1 , for x>0. 

' x=0 1 t=0 1 t=0 


Following the method of multiple scales outlined above, we find it 
convenient to define the differential operators by 

(4.4) = d/dx^ , i = 0, 1, ... . 

We then look for solutions for p and u in the form of (3.2), where, in 

particular, p ^ ^ and u^^ satisfy equations (3.4). Then the general solution 
to equation (3.8) is given by <p = F(t-x) + G(t+x), where F and G are 
arbitrary functions of their arguments and may also depend on the spacial 

scales x^\ x^ 2 \ ... . Since we shall restrict our attention to solutions 
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which propagate only in the positive x direction, we set G = 0 and write 


(4.5) 


, (1 > =u (0 > = fit-,, x (1) , * (2) , ...), 


where the exact form of the function f will be determined from the boundary 
and initial conditions (4.3) of the problem. 

Using the solutions (4.5), equations (3.5) with k=l for and u^ 

become 


(4.6) 


dp {2) Id t + D Q u (1) = -^f + D q ( f 2 ) ) 
du (1) /dt + D Q p (2) = -(d^ + ( ( 7 -I) / 2 )D Q (f 2 )j . 


The general solution to (4.6) for p^ 2 ^ and can be written as 


(4.7) 


p {2) = -xfbjf + ((l+ 7 )/ 4 )D 0 (f 2 )] + f 1 (t-x, x (1) , x (2) , ...), 

u (1) = -x^f + ((l+ 7 )/ 4 )D 0 (f 2 )] 

+ (( 7 - 3 )/ 4 )(f 2 ) + f 1 (t-x, x (1) , x (2) , ...). 


In (4.7), f^ is an arbitrary function of its arguments, which will eventually 
be determined by the initial and boundary conditions of the problem. 

In order for p 1 ' 2 '* and u^ to grow no faster that p ^ and u^°\ 
respectively, as x becomes large, we must require that the term in brackets in 
equations (4.7) vanishes, i.e. 


(4.8) 


D x f + ((l+7)/4)D 0 (r) = 0 


Condition (4.8) is an equation to determine the 


"x (1) behavior" of f and, 
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hence, the x^ behavior of both and Equation (4.8) can be 

expressed in terms of Burgers’ equation, with ((l+'y)/4)x^^ playing the role 
of the "time" variable and x playing the role of the spacial variable. In 
particular, the solution to (4.8) which satisfies the "initial" condition that 

f = f(t-x) when x^=0 can be expressed as 


(4.9) 


f = f(r), where r = t-x + (('y+l)/2)x^f(r) . 


As an application of these results, we let the (dimensional) fluid 
velocity at x=0 be u = Usin(wt). We then use U as the typical velocity of the 
flow and let L = c^/w be the typical length associated with the flow. (Note 

that k = 1/L = w/Cq is the linear wave number for this example.) Then the 

boundary condition for the flow, when expressed in terms of the nondimensional 
variables (2.4), becomes 

(4.10) u = sin(t), at x=0 , for all t^O (boundary condition). 


Using the expansions (3.2) in (4.10) we find that the individual coefficients 
(k) (k) 

u and p satisfy the initial and boundary conditions 


(4.11) 


(4.12) 


u (k >= 0 

/> (k+1) = 0 


u <0) = sint 
u (k) = 0 


at t=0 (for all x>0) for k = 0, 1, 2 


at x=0 (for all t^O) for k = 1, 2, 3, 


From equations (4.5), we find 


(4.13) 


u (0 > . , (1) = m-x, x (1) , *< 2 \ ...I 
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If we were to terminate our perturbation expansion at this point and set x 

for i^l equal to zero, then the conditions (4.11)— (4.12) would yield the 
results 


(4.14) 



sin(t-x), for O^x^t , 
0 , for x>t . 


This, of course, is just the classical wave solution of linear acoustics. 
However, if we continue the perturbation solution outlined above and use the 
solution (4.9), we find 


(4.15a) 


(4.15b) 


u (0) = , (1 > 


sin(r) , for 0 ^x<t , 
0 , for x>t , 


where 


r = t-x + ( ( 7 + 1 ) / 2 )x^ 1 ^sin(r) . 


P 


( 2 ) 


With 
and u 


p (1) 

a) 


determined from (4.15), we find that the expressions (4.7) 
which satisfy conditions (4.11) and (4.12), reduce to 


for 


(4.16) 


p ^ = ((3-qf)/4)sin 2 (t-x) , 

u^ = ( ( 7 - 3 ) /4) ^sin 2 (r) - sin 2 (t-x)j, 


for O^x^t, with p^ 2 ^ = u^ =0 for t>x. 


Combining the expressions we have obtained for p'^ and p^\ we can 


write 


(4.17) 


p = 1 + Msin(r) + M 2 ( ( 3 - 7 )/ 4 )sin 2 (t-x) + 0(M 3 ) + 0(M 4 x) , 
u = sin(r) + M( ( 7 - 3 ) /4) ^sin 2 (r) - sin 2 (t-x)j 
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+ 0(M 2 ) + 0(M 3 x) , 

for O^x^t, and p= 1 and u=0 for x>t. Here t is defined by equation (4.15b). 

To illustrate these results, in Figure 1 we have plotted approximations 
to (p- 1) as a function of distance from the source at a fixed time. The 
multiple scales approximations (solid lines) were obtained using the first 
three terms on the right side of the first equation in (4.17) and are plotted 
for two different values of M. In Fig. 1(a), M=0.01 and the resulting wave 
lies close to linear solution (4.14) (short dashed lines), with only a small 
distortion due to the nonlinearities in the governing equations being visible. 
In Fig. 1(b), M has been increased to 0.04. For this case, the gradual 
steepening of the wave as x increases is evident, with the wave apparently 
approaching an "N-wave" as x increases. An approximation based on Whitham’s 
(1974) first order approximation, obtained using essentially the method of 
strained coordinates, is very close to (but not identical to) the multiple 
scales solution, and hence is not shown on these plots. In these figures we 
have also plotted a solution (circles) to the basic equations (4.1)-(4.2) 
obtained by purely numerical means. This numerical solution was obtained by 
writing the basic equations in characteristic form and then using a MacCormack 
predictor-corrector technique on the resulting equations. (We thank Dr. 

Willie R. Watson of the NASA Langley Research Center for carrying out these 
calculations for us.) As the figures illustrate, there is very good agreement 
between the multiple scales solution and the numerical solution for these 
values of the Mach number. We shall comment further on these results in 
section 9. 

5. Three-Dimensional Flow - The Vibrating Sphere 

We now wish to describe the acoustic field in an inviscid, isentropic 
fluid outside a spherical region, when the normal component of the fluid 
velocity is specified on the surface of this region. It is convenient to 
think of this spherical region in one of two ways. First, we may regard it as 
an actual sphere, whose surface executes prescribed vibrations, with the 
boundary conditions applied on the time averaged surface of the sphere. Under 
this interpretation, the results which follow are extensions of several 
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"classical" problems in acoustics. Alternatively, we may regard the region as 
a (mathematical) sphere which encloses a bounded region of "complicated" fluid 
flow. With this interpretation, our results may be applicable to the problem 
of determining the acoustic field radiated from bounded (perhaps turbulent) 
fluid flows. In particular, our results will show that only the normal 
component of the flow velocity needs to be specified on this spherical surface 
in order to determine the radiated acoustical field. 

To begin, we let the origin of a (nondimensional ) Cartesian coordinate 
system (x, ,x x ) coincide with the average position of the center of the 

sphere, whose average radius is a constant a. We then introduce spherical 
coordinates (r ,0,t) (see Figure 2(a)) related to (x^^) by the relations 

x = rsin(0)cos(i/O , x 0 = rsin(0) sin(i/0 , x = rcos(0). (Here all lengths have 

been nondimensionalized by referring them to the average radius a of the 
spherical region.) Then, in terms of these coordinates, equations (2.5)-(2.6) 
become 


(5.1) dp/d t + M((l/r 2 )(d/dr)(r 2 pu r ) + (l/rsin(0) ) (d/dd) (sin(^)pu 0 ) 

+ (l/rsin(tf))(9/^')(pu^ r )j = 0, 

(5.2) dp/dx + p 2 ^Mdu^/dt + M 2 ^u^Ou^/Sr) + (u^/r) {dn^/d0) 

+ (u^/rsin(0)) {da^/d^) - (u^) 2 /r - (u^) /rjj = 0, 

(5.3) dp/dO + rp 2 ^McJu^/dt + M 2 ^u^. ( dn^ /dr ) + (u^/rMdu^/d#) 

+ (u^/rsin(0) ) (du^/ety) + (u r u^)/r - (u^) cot(0)/rjj = 0, 

(5.4) dp/di + (rsin(0))p 2 " 7 |M&ydt + U 2 {xi^du^/ dx) + (u^/r) {dn^/dO) 


In (5.1)— (5.4) , u r , u^, represent the fluid velocity components in the 

positive r, 0 and ^ directions, respectively. 

We shall assume that u r is a specified function, say V(0,^,t), at r=l for 

all t^O , and that u^, u^, u^ and p-1 are all specified to be zero at t=0 for 
all l<r<M, i.e. 


(5.5) 


u 


r 


= v(fl,vM), 

r=l 

t^O 


( WV 


t=0 

r>l 


(0,0,0) and 


t=0 

r>l 


1 . 


Following the method of multiple scales, we introduce the spacial scales 
r., i = 0, 1, 2, ... , related to r by 


(5.6) 


= M r, i = 1, 2, 


Then, for small values of M, we look for solutions for the density and 
velocity components in the form 


(5.7) 


co 


p = 1 + M/> (1) + M V 2) + ... = Y2 

j=0 


, with p^^= 1, 


00 


u r = u<°> * Mu<‘> * M 2 uJ 2 > + ... = 

j=0 


with analogous expressions holding for u^ and u^. Here each of the 
coefficient functions p ^ ^ , etc., is independent of M, but, in general, 
will depend upon t and the spacial scales r., i.e. = 

/> (j) (t,r 0 ,r 1 ,r 2 ,...), } = u< J } ( t , r Q , ^ . . . ) , ... . 

To determine these coefficient functions, we substitute the expansions 
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(5.7) into equations (5. 1)— (5.4) and use the relation 

(5.8) did r = D Q + W> 1 + M 2 D 2 + ... , where Dj = d/dr. . 

We then collect coefficients of like powers of M on the left side of each 
equation, as described in section 3, and find the following system of 

equations satisfied by uj.^ , u^^ , an ^ : 

(5.9a) dp (1) /dt + (l/r 2 )D 0 (rV 0) ) + ( 1 /rsin(0) ) (d/d6) ( sin(0 )u^ 0) ) 

+ (l/rsin(0) )du^ 0 ^ ld$ - 0, 


(5.9b) 

D 0 p (1) + du[ 0 ) /5t = 0, 

(5.9c) 

dp^^ IdO + rdug®^ Idt = 0, 

(5.9d) 

dp (1) /ty + rsin(0)du^ O) /dt = 0 ; 

(5.10a) 

<9p (k+1) /dt + (l/r 2 )D 0 [r 2 uJ k) ] 


+ (l/rsin( 0 ))(d/d 0 )(sin( 0 )u 0 k ^) + (l/rsin( 0 ))du^ k )/cty 

(5.10b) 

(k+ 1 ) a(k) /;) , 

D q p + du f Id t = G r , 

( 5 . 10 c) 

3p (k+1) Id6 + rdM^/dt = G^ k) , 

(5. lOd) 

dp ik+1) /ty + rsin( 0 )du^ k) Idt = G^ k) , for k=l, 2 


(Equations (5.9) follow from the terms in equations (5.1)-(5.4) which are 
O(M), while equations (5.10) follow from the terms in these equations which 

are 0(M k+1 ).) Here the functions F (k) , G^ k) , G^ k) and G^ k) depend only upon 
with j<k+l and u^\ u^ ) , and with j<k. 

In the next sections we shall show how equations (5.9)-(5.10) can be 
solved recursively. 
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6. Solution for the lowest order perturbation coefficients 


The form of equations (5.9) suggests that we look for a solution for the 
lowest order perturbation coefficient functions in the form 

(6-D = -d<p/d t, u£^ = dp/dr, 

u< 0) = (l/r)dp/dO, uj° } = (l/rsin (0))dp/fy, 

where p is a smooth function of its arguments, which still needs to be 
determined. Using equations (6.1), we see that equations (5.9b)-(5.9d) are 
satisfied for any choice of p, while equation (5.9a) leads to the requirement 
that p must satisfy 

(6.2) d 2 p/d t 2 = (l/r 2 )(d/dr) {r 2 dp/dr} 

+ (l/r 2 sin(0))(d/d0) ^sin(0)dp/dtfj + (l/r 2 sin 2 (0))d 2 p/d^ 2 , 

which is just the usual (linear) wave equation. (Here we have used the 
relation D Q a d/dr Q = d/dr.) This equation has solutions of the form 


(6.3) p = F n (r, t)P^(cos(0))cos(m^-/?) , for n = 0, 1, 2, ... , 

and m=0, 1, ... , n, 

where is the associated Legendre polynomial, (3 is an arbitrary constant, 
and satisfies the differential equation 

d 2 ? Jdt 2 = (l/r 2 )(d/dr)(r 2 <9F /dr) - (l/r 2 )n(n+l)F . 

11 k n / n 

The outgoing solution for F can be expressed as 

n 
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(6.4) 


F (r, t) = y ' a ,f (n+1 j) (t-r+l)/r^, 
n L . n,j 

3=1 

>n,l “ - 1 ’ a n, j + 1 - (("+H)U*n)/<2j))a nJ , j = 1, 2 n, 


where f is an arbitrary (smooth) function of its argument. 
If we now apply the boundary condition (5.5) 


V(0,^,t) = u f (r ,6,^, 


‘>| = ^M j u‘ j) (r 0 , r r 

r - 1 j=0 


• • t) 


we see that we can set 


(6.5) 


u; w/ = V(6,t, t), for j=0, 
r=l 


= 0 , for j^l . 


We shall assume that V can be expressed as (a linear combination of terms of 
the form) 


( 6 . 6 ) 


V(0,^,t) = Q(t)P™(cos(0))cos(mi/ r -/?) , 


where Q(t) is a specified function of t, and n and m are non-negative 
integers, with 0^m<;n. Using this expression for V, along with the condition 

(6.5) on u^\ as well as the expressions (6.3) for <p and (6.1) for we 

find that f must satisfy the condition 


(6.7) f (n+1) (t) - Q(j(j+l)+n(n+l))/(2j))a n j f (n+1 j) (t) 

j=l 
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- <n+1)a n,n + l fU) ' « (t) ' 


Equation (6.7) is an inhomogeneous, linear ordinary differential equation of 
order n+1 , with constant coefficients, for the unknown function f(t). For a 
few small values of n this equation becomes 


n=0 

n=l 

(6.8) n=2 

n=3 
n=4 


f’+ f = Q(t); 

f”+ 2f ’ + 2f = Q(t); 

f’” + 4f”+ 9f’+ 9f = Q(t); 

f””+ 7f’”+ 27f”+60f’ + 60f = Q(t); 

f””’ + llf””+ 65f”’+ 240f”+ 525f ’+ 525f = Q(t). 


Once f has been determined, we can express the solution for <p as 


(6.9) 


<P = 


' n+1 

r'a . f (n+ l-j)(f)/ r j 

L i n , j 

j=l 


P^(cos(0))cos(m^-/?) , 


where t’ = t-r+1 , In general, f will also depend on the "slow" spacial scales 
Tj , Tg, although this dependence has not been denoted explicitly. 

As an illustration of these results, we consider the special case when 
Q(t) = sin(wt). Then the steady state (periodic) solution to (6.7) (which 
can be viewed mathematically as a particular solution to this equation) can be 
expressed as 


(6.10) f(t) = Asin(wt-a), when Q(t) = sin(wt), 

where, for a few small values of n, A and a are given by 

n=0: A = (l+w^) a = tan *(w), 

n=l : A = (4+u/*) , a = tan * j^2w/ (2-w^) j , 
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n=2: A = (81+9w 2 -2u 4 +w 8 ) ^ 2 , a. = tan * |w(9-w 2 ) /(9-< 4w 2 )j , 

n=3: A = (3600+360w 2 +9w 4 -5w 6 +a; 8 ) _1/2 , 

a = tan -1 (w(60-7w 2 )/(c/-27w 2 +60)) , 

9 A fi ft 1 0 —1 / 2 

n=4: A = (275625+23625W +900w -5w -9w +w ) 

a = tan _1 ^(525-65w 2 +a; 4 )/(525-240w 2 +llw 4 )j . 

Solutions to the homogeneous version of equation (6.7) are of the form f(t) = 

e^ 4 , where the real part of A is negative. Hence the corresponding solutions 
for f(t) are not periodic in t and, in fact, they all decay to zero as time 
increases . 

Before we consider the higher order perturbation coefficients, we can use 
the results above to calculate a first approximation to the radial component I 
of the sound intensity, defined by 

(6.11) I = p 0 c 2 U<Mp*u r >, 

where the symbol "<g>" denotes the time averaged value of the quantity g. 

Using the expressions (6.1), (6.3) and (6.4) above we find, for the case when 
Q is given as in (6.10), that 


I = PgCQU 2 w 2m+2 A 2 |p^(cos(0))cos(m^-/?)j /2r 2 + 0(l/r 3 ), as r-xx. 


For w s ka « 1, we find that the maximum intensity of I, which we denote by 
I , for different values of n, is given by 


n=0: I Q = [/ , o c {)U 2 (ka) 2 /2j (a/r) 2 + 0((ka) 4 /r 2 ) + 0(l/r 3 ); 

n=l: 1^ = |pQCQU 2 (ka) 2 /2j (a/r ) 2 ^(ka) 2 /2 2 j + 0((ka) 8 /r 2 ) + 0(l/r 3 ); 
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n=2 : I 2 = ^ 0 c 0 U 2 (ka) 2 /2j(a/?) 2 (^(ka) 4 /3 2 ] +0((ka) 8 /r 2 ) +0(l/r 3 ). 


In general, for any positive integer n, the ratio of the maximum intensity I 

~ 2n 2 

to the maximum intensity Iq of the pulsating sphere approaches (ka) /(n+1) , 

as ka approaches zero. 

7 . Solutions for the higher order perturbation coefficients 
Equations (5.10) with k=l can be written as 


(7.1a) 


(7.1b) 


(7.1c) 


(7. Id) 


dp (2) /at + (l/r 2 )D 0 (rV 1) ) 

+ (l/rsin(0)) {dldO) (sin(^)u^ 1 } ) + (l/rsin(0) )du^ } /dvf- = 
-[(l/r 2 )D 0 (r 2 /) (1) u[ 0) ) + D lU ^ 0) 

+ (l/rsin(0)) {d!d6) (sin(d)p^^u^^ ) 

+ (l/rsin(0))(5/^)(/j (1) u| O) )), 

d n p (2) + du (1) /at = 

0^ r 

-^ 0) D 0 u[ 0) + (2-7)p (1 W 0) /dt + D x p (1) 

+ (u^ O) /r)(du|; O) /d0) + (u^ O) /rsin(0))(du[ O) /^) 

, (0k2, , (0) , 2 , \ 

- {u 0 ) / r - (u^ ) / rj, 

dp^ 2 ^ /dO + rdu^^/dt = 

-r((2- 7 )p (1) ^ 0) /at + u[ 0) D 0 4 0) 

+ (u^ O) /r)(0uJ O) /d0) + (u^ 0) /rsin(^) ) (du^ 03 /cty) 

+ (u r )/r - (u^ ) cot(0)/rj, 

dp i2) /d+ + rsin(0)5u^ 1) /5t = 
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-(rsin(0) ) \{2-q) p^ I dt + u^D^ 

+ (u^/r) (du^ 0 ^ IdO) + (u^ O Vrsin(0))(5u^°^/d^) 

, (0) (Oh, , (0) (0), ., a . , 'i 

+ (u r )/r + (u^ u^ )cot(0)/rj. 

To begin our construction of solutions to equations (7.1), we first 
examine the behavior of the right sides of these equations as r-»oo. Using 
equations (6.1), (6.3), and (6.4) we find 

= -dp/d t = (l/r)f ^ n+ ^P™(cos(0) )cos(m^-/?) + 0(l/r 2 ), 

(7.2) u£°) = dptdx = (l/r)f^ n+1 ^P™(cos(0))cos(i#-/?) +0(l/r 2 ), 

uj 0) = {l/r)dp/d6 = (l/r 2 )f ^P^’ ( cos(0) ) (sin(0) )cos(i#-/?) +0(l/r 3 ), 
u< 0) = (l/rsin(0) )d<p/d4' 

= (m/r 2 sin(fl) )f ^ n ^P™(cos(0) )sin(m^-y3) +0(l/r 3 ), as r-*x. 

(Here it is understood that the argument of f is t’=t-r+l.) Then, using these 
expressions, we can write equations (7.1) as 

(7.3a) dp i2) /d t + (l/r 2 )D 0 (rV 1) ) 

+ (l/rsin(0 ) ) (d/d&) ( sin(0)u^ 1 ^ ) + (l/rsin(0) )du^ 1 ^ /cty = 

-DjU^ + (2/r 2 )f^ n+1 ^f^ n+2 ^F^(cos(tf))cos(m^-/?)j +0(l/r 3 ), 

(7.3b) D q p ( 2) + duj^/dt = 

2 

+ ((7-l)/r 2 )f^ n+1 ^f^ n+2 ^P™(cos(fl))cos(mi/'-/3)j +0(l/r 3 ), 

(7.3c) dp^ 2 ^ /dO + rdai^/dt = 
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( (I- 7 ) /r 2 ) P®(cos(0))P^’(cos(0))sin(0)cos 2 (m^-/?) + 0(l/r 3 ), 

(7.3d) dp^ !di> + rsin(0)du^^ /dt = 

2 

m( ( 7 -I ) /r 2 ) |f ( n+ l)p®( cos( 0 ) ) j cos(i#-/?)sin(m^-/?) + 0 (l/r 3 ). 

We now examine the effect on the solutions p^ 2 \ u£^, etc., of each of 
the terms on the right sides of equations (7.3). In particular, we must first 
examine the behavior of the terms DjP^ and Djuj; 0 ^ as r^oo. Now, if we assume 

that D^p^ and D^u^^ are 0(l/r) as r-*x>, then these terms will give rise to 

terms in p^ and u^ which are 0(log(r)) as r->oo. However, since p^ must 

decay at least as fast as p^ , which is 0(l/r) as r-x», we see that this 
cannot be allowed. Consequently, we shall now require that D^p^^ and D^u^^ 

2 

are 0(l/r ) as r-»co, and we will show that this assumption allows us to 
determine the quantities of interest in a consistent manner. Using (7.2), we 
can express this assumption as 

D l p(l) = D l u r° )+ 0(1/r3) 

= (1/r 2 ) (rDjf ^ n+1 )jp^(cos(0))cos(m^-/?) +0(l/r 3 ), 

where rDjf^ n+ ^^ is bounded as r-xx. Also, as we shall show below, it is 

convenient to think of the argument t’ of f as being replaced by a new 
argument, say, t = r( t-r+1 , 6 ,i, ^ , r 2 , . . . ,M) , which must be determined. Then 

we can write 

(7.4a) D 1 f (n+l) (r) = f (n+2) (r) (D^) . 
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Furthermore, the angular dependence of the terms on the right sides of 
equations (7.3) suggests that should be proportional to p“(cos(0))« 
cos(m i-fi) . Thus, we write 


(7.4b) 


D , T = h(r)P m (cos(0) )cos(m^-/3) , 
1 n 


where h(r) must be determined. With these assumptions, we see that all of the 
terms on the right side of (7.3) are 0(l/r ) (at least) as r-»oo, as required. 
Thus, we can express the (particular) solutions for \ \ , and u^ 

in the form 


.( 2 ) 


u 


.( 1 ) 


(7.5) 


u 


( 1 ) 


u 


( 1 ) 


|(l/r) log(r)Fp - (1/r 2 ) log(r)Fp + (l/r 2 )F 1 | 

‘ (p^( co s(0) c °s(m^-/?)j + 0(l/r 3 ), 

|(l/r)log(r)F^ + (l/r 2 )F 2 | 

' (Vj(cos(0)cos(m*-/?)} 2 + 0(l/r 3 ), 
|(l/r 2 )log(r)F 0 + (l/r 3 )F 3 j 

'2P ra (cos(0))P m ’(cos(0))sin(0)cos 2 (m^-/3) + 0(l/r 4 ), 

n n 

j(l/r 2 sin(0))log(r)F o + (l/r 3 sin(0) )F 4 j 

f \ 2 4 

’2mfp^(cos(0) )j cos(m^-/?)sin(m^-/3) + 0(l/r ), 


where each F. = F.(r,r) is a function to be determined and the primes denote 
J J 

differentiation with respect to r. Substituting the expressions (7.5) into 
equations (7.3), and using (7.4), we obtain the relations: 


21 


(7.6a) Fq = (( 7 +l)/ 2 )f (n+ 1 ) f (n+2) - rf (n+ 2 ) h(r) + 0(l/r), 

(7.6b) F’ = F’ + (( 3 - 7 )/ 2 )f (n+ 1 ) f (n+2) + 0(l/r), 

2 

(7.6c) F’ = F 2 + ((7-37)/2)[f (n+2) ] + 0(l/r), 

2 

(7 . 6 d) F’ = F 2 + (( 7 - 7 )/ 2 )(f (n+1) j + 0(l/r). 

In (7,6), the function can be determined by examining the terms in 

3 

equations (7.3) that are 0(l/r ). From the form of the solution (7.5), we see 
that we must demand that F Q = 0 (l/log(r)), as r-*oo, since otherwise p ^ will 

not decay as fast as the 0(l/r) rate of decay of p^^ . Using this condition, 
we see that the right side of equation ( 7 . 6 a) must be 0 (l/log(r)), as r-*». 
Using the relation (7.4b), this condition can be written as 

(7.7) D x (g) + ((7+l)/4)pJJ(cos(^))cos(m^-/?)(l/r 2 )D 0 (r 2 g 2 ) = 0(l/(r 2 log(r)) ) , 

as r-*oo, where g s f^ n+ 1 ^/r. (Here we have assumed that D^r = -1 + 

0(l/log(r)) and dr/dt = 1 + 0(l/log(r)), as r->oo. See equations (7.9) below.) 
Equation (7.7) is the major result of this section and it serves to 

determine the behavior of p^\ u£°\ u^, and u^ as functions of ^ . It 
can be interpreted as a "spherical" Burgers’ equation in the "time-like" 
variable ( ( 7 + 1 ) /2)P^( cos( d ) ) cos(mi^-/3)r^ . In particular, we can express the 

solution of (7.7) which satisfies the "initial" condition g = f ^ n+ 1 ^ ( t+l-r ) /r 
when M=0 as 

(7.8a) g = f^ n+ D(r)/r, where 
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(7.8b) 


t = t-r +1 + ((7+l)/2)F n (r,r 1 )P™(cos(0))cos(m^-/3) » 


n+1 r 

F (r , r 1 ) = -Y " a . f (n+2 ~j } (r^ f 1 (l/s) j ds 
n 1 L— i n,j J M 

j=l 


n+1 


= Mf 


(n+1) (r)log(r 1 /M) + a n j f (n+2 j } (r) ((M/^ ) j 1 -l) 

j=2 


where a . = a - / ( j — 1 ) , j = 2, 3, ... , n+1. 
^ » J “ J J 


Thus, both the and behavior of p^\ uj. \ \ and are completely 

determined by equations (6.9), (6.1), with t’ replaced by r, and (7.8b). 

Before we consider some specific applications of these results, we note 
that, from its definition in (7.8b), r has the following properties: 


(7.9a) t •* t-r+1, as M+0; 

(7.9b) t -» t-r+1, as r-»l; 

(7.9c) dr/dt = 1 + 0(l/log(r)) , dr/d r = -1 + O(l/log(r)), as r^oo; 
(7.9d) dr/dO = O(M) , and dr/d ^ = O(M) , as M+O; 


and 


(7.9e) 


D, r = ((7+l)/2)P m (cos(0))cos(m\M)f (n+1) (r)/r + 0(l/r 2 ), as r-*x. 
1 n 


Properties (7.9a) and (7.9b) insure that the linear solution is recovered, 
either as M+0 or as we approach the surface of the sphere. Property (7.9c) 
shows that our definition of t is consistent with the assumptions made above, 
while property (7.9d) shows that the angular variation of r is "small" and, 
hence, will be described by the next level of the perturbation expansion. 
Property (7.9e) shows that is 0(l/r), as r-*x>, and hence t satisfies the 

requirement that rD^r is bounded as r-»oo. We shall comment further on the 

properties of t in section 9. 
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8. Appl i cations 


a) The Pulsating Sphere 

As our first application, we consider the pulsating sphere (see Figure 
2(b)). In this case, the (nondimensional ) radial component of the velocity of 
the surface of the sphere is given by V = sin(wt) and hence we set n=0 and m=0 
in our results from sections 6 and 7. Then the steady state solution for f(t) 
is given by (6.10) with n=0, and ip is given by (6.9) with n=0, i.e. 

(8.1) f(t) = (1+w 2 ) ^sin(wt) - wcos(wt)j and <p = -f(r)/r, 

where r is determined by (7.8b) with n=m=0, i.e.. 


(8.2) t = t-r+1 + ((7+l)/2)Mf’(r)log(r 1 /M). 

Using these definitions in equations (6.1) and (7.1) we find for this case 
p^ = |w/ ( (l+c</*)r)j (cos(o;r) + wsin(wr) ) , 

( 2 ) 

p = 2bgWcos(2wr) /r - 2b 2 wsin(2wr) /r 

- l/(4r 4 (l+w 2 )) + (l-7)w 2 /(4r 2 (l+w 2 )) 

+ ((17+7)w 2 /(16r 3 (l+w 2 ) 2 ) - ( 13-37) (w 2 -l)/(64r 4 (l-HJ 2 ) 2 ) 

(8.3) + (5+7)w 2 (u 2 -l) /(8r 2 (l+w 2 ) 2 ) j cos(2a;r) 

+ [(13-37)w/(32r 4 (l+« 2 ) 2 ) - (5+7)w 3 /(4r 2 (l+a; 2 ) 2 ) 


+ (17+7)w(w 2 -l)/(32r 3 (l-H*; 2 ) 2 )j sin(2a;r) 

+ 0(l/r 5 ), 

and 

b 2 = -3 w 2 (14-27+17w 2 +tw 2 )/(16(1+ci> 2 ) 2 (1+4w 2 )) , 
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b„ = w(29-37-17w 2 +15'7W 2 -64w 4 )/(32(1+w 2 ) 2 (1+4w 2 )) . 
o 

(The corresponding expressions for the perturbation coefficient functions 

and u 4 for the radial component of the velocity are listed in the 
r r 

appendix . ) 

In Figure 3, we have plotted r(p^^ + Mp^ 2 ^) (solid lines) as a function 
of r at a fixed time for w=1.25 and for two different values of M. Here it is 

convenient to think of p^^ + Mp^ 2 ^ as our approximation to the normalized, 
nondimensional acoustical density field (p-l)/M (from equation (3.2)). In 
Figure 3 we have also plotted the approximation to r(p-l)/M based on Whitham’s 
(1974) first order solution (long dashed lines), the classical linear solution 
(short dashed lines), and a solution obtained by purely numerical methods 
(circles). This numerical solution was obtained using a MacCormack 
predictor-corrector scheme on equations (5. 1)— (5.4) . As Figure 3 illustrates, 
there is very good agreement between the multiple scales perturbation 
solutions and the numerical solutions, even for a Mach number of 0.3, where 
the gradual steepening and asymmetry of the wave is apparent. 

b) The Oscillating Sphere 

As a second application, we consider the oscillating sphere (see Figure 
2(c).) In this case, the radial component of the velocity of the surface of 
the sphere is given by V = cos(0)sin(wt) and, hence, we set n=l and m=0 in our 
results from sections 6 and 7. Then the steady state solution for f(t) is 
given by (6.10) with n=l , and ip is given by (6.9) with n=l, i.e. 

(8.4) f(t) = (4+w 4 ) -1 ^(2-w 2 ) sin(wt) - 2ucos(wt)j , 

<p = - cos(0) ff ’ (r) /r + f(r)/r 2 j, 

where r is now determined by (7.8b) with n=l and m=0, i.e. 
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(8.5) r = t-r+1 + ((7+l)/2)M|f”(r)log(r 1 /M) - f ’ (r) ((M/r^-lj jcos(0) . 

Using these expressions in equations (6.1) and (7.1), and using Wolfram’s 
(1991) symbolic computation system Mathematica, we find for this case 

|w(2-w 2 +2ru; 2 )cos(wr) /r 2 + w 2 (2-2r+rw 2 )sin(wr) /r 2 j cos(0) / (4+w 4 ) , 

p^ = |2w(f 0s cos(2w7-)-f 0c sin(2wr) )/r - (3w/2) (f 2g cos(2wr)+f 2c sin(2w7-) ) /r 3 

+ 3w 2 (f 2(; cos(2a;r)+f 2g sin(2a;r))/r 2 + 2w 3 ( f 2g cos(2wr)-f 2c sin(2wr ) ) /r 

+ (l-7)w 2 (4+w 4 )/8r 4 + (l-7)w 4 (4+w 4 )/8r 2 
+ ^(323+3-y)w 2 (2-w 2 )/64r 5 - (7+7 )c/(2-w 2 ) /2r 3 

- (259+197)w 2 (4-8w 2 +ti; 4 ) /192r 4 + (5+7)w 4 (4-8w 2 +w 4 ) /16r 2 j cos (2wr) 

+ ^( 5+7) w 5 ( 2-w 2 ) /4r 2 - (259+197)w 3 (2-w 2 )/48r 4 

+ (7+7)w 3 (4-8w 2 +w 4 ) /8r 3 

- w(323+37) (4-8w 2 +w 4 ) /256r 5 j sin(2wr)|/(4+w 4 ) 2 
+ cos(20) j(9w/2) (-f 2g cos(2wr)+f 2c sin(2wr) ) /r 3 

- 9w 2 ( f 2c cos(2wr)+f 2g sin(2wr) ) /r 2 + 6w 3 (f 2s cos(2wT)-f 2c sin(2wr) ) /r 

+ (3-7)w 2 (4+w 4 ) /8r 4 + ( (l-7)u 4 (4+w 4 )/8r 2 
+ ^3(33+7 )cj 2 (2-w 2 ) /32r^ - (17+7)w 2 (4-8w 2 +w 4 ) /16r 4 

- 3(9+7)w 4 (2-ci; 2 ) /8r 3 + (5+7)w 4 (4-8w 2 +w 4 ) /16r 2 j cos(2wr ) 

+ ^(5+7)w^(2-w 2 ) /4r 2 - ( 17+7)w 3 (2-a; 2 ) /4r 4 

+ 3(9+7)w 3 (4-8w 2 +w 4 ) /32r 3 + 
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- 3w(33+7)(4-8cJ 2 +a; 4 )/128r 5 jsin(2wr)j/(4-H l ; 4 ) 2 + 0(l/r 5 ). 

( 2 ) 

The constants f , f Q , f 2g , and f 2<j which appear in the definition of p 

are related to the homogeneous solution to equations (7.1) and are defined in 

the Appendix. (The form of the solutions for u^. , u^ , u^ , and are 

also presented in the Appendix.) 

The form of these expressions allows us to make an interesting 
observation concerning the angular dependence of the sound radiated from an 
oscillating sphere. Specifically, it is interesting to examine the density 

and velocity in the 0=7r/2 direction. In this direction, both p ^ ^ and 

are zero , while uj^^ is nonzero. In contrast , for this direction both p 

and uj 1} are nonzero, while is zero (since it contains a factor of 

sin(20) ) . Thus, there is a "small" (0(M 2 )) amount of sound radiated in this 
direction, which has a characteristic frequency of 2w. In particular, in the 
far field we can use the expressions above to write 

p (2) =u (1) + 0(l/r 2 ) 

r 

= (2w/r) [(2f 2s w 2 -f 0g )cos(2wr) + (f 0c -2f 2c w 2 )sin(2ur)j +0(l/r 2 ), 

along 0= 7r/2. 

In addition, we find that the expressions for uj/^ and correspond to a 

rot at ional flow, even though uj,^ and correspond to an i r rot at ional flow. 

The presence of vorticity in the flow can be interpreted as being due to 
tangential acceleration of the fluid at the boundary, as discussed by Morton 
(1984). This vorticity appears only very near the sphere and is given by 
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^xu = i^Msin(20) |j^( 7 +l )w 4 / (4(4+w 4 ) )j ( log(r ) ) /r 3 

+ |w 4 (4-8u 2 +w 4 ) /(4+w 4 ) 2 j cos(2wr) /r 3 
+ ^4w 5 (2-w 2 )/(4+w 4 ) 2 jsin(2a;r)/r 3 + 0(l/r 4 )|. 

In Figure 4, we have plotted approximations to r(p-l)/(Mcos(0) ) as a 

2 

function of r at a fixed time for 0=0 and 0= tt/3, as well as r(p-l)/M at 
0= 7t/2, with w=1.25 and for two different values of M. These approximations 

are the multiple scales solution r (// 1 ^+Mp^ 2 ^ ) /cos(0) (solid lines), the 
classical linear solution (short dashed lines), and a solution obtained by 
purely numerical means (circles). For this case, we found that numerical 
approximations based on a Jameson type Runge-Kutta finite volume scheme, 
which is second order accurate in both space and time, gave better agreement 
with the multiple scales solution than numerical approximations based on the 
MacCormack scheme, as in the case of the pulsating sphere. (We thank Mr. 

David Lockard and Dr. Kenneth Brentner of the NASA Langley Research Center for 
carrying out these calculations for us.) Again, the agreement between the 
multiple scales solutions and the numerical solutions is qualitatively good. 

c) A "Squishing" Sphere 

As a third example, we consider the case when the radial component of the 

2 

surface velocity of the sphere is given by V = sin (0)sin(2^)sin(tt;t) (see 
Figure 2(d)), Thus, we set n=2 and m=2 in our results of sections 6 and 7, 
and find from (6.10) with n=2 that the steady solution for f(t) is given by 

(8.7) f(t) = (3(81+9 w 2 -2w 4 +w 6 ) ) _1 |^(9-4w 2 ) sin(wt) - w(9-u> 2 )cos(wt)j , 
while the solution for p is given by (6.9) with n=2, i.e., 

(8.8) <p = -^f”(r)/r + 3f’(r)/r 2 + 3f(r)/r 3 j3sin 2 (0)sin(2^) , 
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with t determined by (7.8b) with n=2, i.e., 


t = t-r +1 + ( ( 7 + 1 ) / 2 )F 2 (r ,r 1 ) 3 sin (0) )sin( 2 f) , 
F 2 (r, ri ) = Mjf’”(r)log(r 1 /M) - 3 f ” (r) [(M/r^-lj 


- (3/2)f’(r)[(M/ ri ) 2 -l] 


Using these expressions we find 
p^ 4 ^ = sin( 2 ^) sin 2 ( 0 ) ^ 

( (27w-12w 3 ) / r 3 +(27w 3 -3u/*) /r 2 -(9w 3 -4u/*) / r )cos(wt ) + 

((27w 2 -3w 4 )/r 3 -(27w 2 -12w 4 )/r 2 -(9w 4 -w 6 )/r)sin(a;T)j / (81+9w 2 -2u 4 +w 6 ) , 

while expressions for and u^^ , as we ll as the form of the next 

order perturbation coefficient functions, are given in the Appendix. 

2 

In Figure 5 we have plotted approximations to r(p-l)/(Msin (0)sin(2^)) 
as a function of r at a fixed value of time with M= 0.6 for ^ = n/ 12 , n/6, and 
■kI 4 at (a) 0=7r / 4 and (b) 0=tt/2, using the multiple scales solution 

rp /(sin ( 0 )sin( 2 ^)) (solid lines) and the classical linear solution 
(dashed lines). Due to the excessive storage requirements for this fully 
three-dimensional, time dependent problem, no numerical approximations were 
computed. 

9 . Conclusions and Discussion 

In this section we shall first make some general observations concerning 
the perturbation approach we have employed to the problem of sound generation 
by vibrating bodies, and then comment specifically on the insights gained from 
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the examples we have presented. 

We note first that the multiple scales perturbation approach has allowed 
us to "capture" many of the salient nonlinear features of acoustic wave 
motion. For example, as the Mach number increases, the nonlinear 
characteristics of the gradual steepening of the waves, the asymmetries due to 
"second harmonic” terms, as well as the convergence of the wave to an "N-wave" 
profile, are all evident in the examples we have considered. These phenomena, 
of course, are not predicted by linear theory. Thus, it appears that the 
method has allowed us to obtain approximate analytical solutions which are 
valid for (much) larger values of M than are the classical linear acoustic 
solutions. In addition, the solutions assume a rather simple form, being 
essentially of the same form as the classical linear solutions, but with a 
different argument. That is, the classical retarded time t+l-r in the linear 
solution has simply been replaced by r, where r is defined implicitly. Other 
perturbation approaches to this problem have resulted in more involved 
expressions for the quantities of physical interest (cp., e.g., Crow (1970)). 
(We should note that, if the present method of analysis were to be continued 
to a higher order than that presented here, various terms of the form 



should be added to the expansions (5.7). 


However, since the 


2 

present analysis was terminated at terms which are 0(M ), these terms were not 
needed and, hence, were not shown explicitly.) 

The general form of the equation for r also provides an insight into the 
interplay between the effect of the Mach number and the angular variation of 
the acoustical source on the form of the acoustic wave. For example, if M=0 

or if the term P^j(cos(0) )cos(m^-/3) is zero (corresponding, say, to a 

particular value of 6 or ^) , then from (7.8b) we see that r- t-r+1 and, hence, 
our solution (to leading order) reduces just to the linear acoustic solution. 


However, if M^O, but the angular term P^( cos (0) ) cos (i#-/9) is zero , then our 

solution (to leading order) again reduces to the linear solution and hence the 

2 

nonzero Mach number has only a second order (i.e. 0(M )) effect on the form of 

the solution. These observations motivate us to define an "effective" Mach 

number M by 
et t . 
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M g ^ = MP™(cos(0))cos(m^-/?) . 

Thus, it is the (angularly dependent) effective Mach number M e ^ , and not 

just M, which determines the magnitude of departure of the wave from its 
linear form. This phenomena is illustrated in Figure 4 for the example of the 
oscillating sphere, for which n=l and m =0 in the formulas above, and also in 
Figure 5, corresponding to n =2 and m= 2 . In Figure 4(a), for example, the 
linear form of the wave is clearly seen, even for nonzero Mach numbers, when 
cos( 0 ) approaches a value such that M^^, is zero. 

In comparing our results with Whitham’s (1974) results, we note that for 
spherically symmetric disturbances (for which n=0) our equations (7.8) reduce 
to his results. However, in more general problems (for which n>0), our 
results differ from Whitham’s equations. In the far field, where the log term 

in the expression for F^ dominates the remaining terms, equation ( 7 . 8 b) 

reduces to 

r = t-r +1 + ( ( 7 + 1 ) 12) f (7-) log(r j/M)MP™(cos(0) )cos ( 1 #-/?) . 

For n=0, this equation is equivalent to Whitham’s equation for r. For n>0, 
we use r^/M = r to note that this equation is the same as Whitham’s result 

when we replace M in his formula by M g ^ . 

For n>0, the definition of t involves terms, through the definition of 

F , which are negligibly small in the far field. However, these terms 
appear to have an interesting (and vital) interpretation. In particular, the 
terms other than the log term in F n give rise to terms in D^r which are 
2 

0 (l/r ) (at least) as r-*», which, in turn, give rise to terms on the right 

3 

sides of equations (7.1a,b) which are 0(l/r ), as r-*». However, without these 
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terms present, no steady state, periodic ( particular ) solution exists for 

p^\ and ul^. To see this, we note first that, if a periodic 

r p y 

solution with period, say, T, for these quantities does exist, we can 
integrate equations (7 . lb)— (7 . Id) over a complete period in t and, with a 
little manipulation, obtain the relations 


T 

0 = (d/d0)J fright side of 
0 
T 

0 = (d/etylj fright side of 


T 

(7,lb)jdt - DqJ fright side of 
0 
T 

(7 . lb) jdt - DqJ fright side of 
0 


(7 . lc ) j dt , 
(7 . Id) j dt . 


The first term on the right side of each of these relations involves D^p , 

while the second terms in each relation is independent of D^p^ \ For 

example, in the case of the oscillating sphere, the first of these relations 
yields the requirement 


0 = (a 0 + l)(7+l)w 4 /(8r 3 (4+w 4 )) . 

Thus, if a periodic solution is to exist, we must have 2 = 4 = a l 2’ as 

stated in equation (7.8b). Similar equations hold for n£2, which led to the 

particular definition of F n which appears in equation (7.8b). 

It is also interesting to note how the distortion of the acoustical wave 
from a linear wave varies with the parameter n. In Figure 6 we have plotted 

rp (1) » r(p-l) /M at M=0.3 for (a) the pulsating sphere (n=0), (b) the 
oscillating sphere at 0=0 (n=l), and (c) the squishing sphere at 0=ir/2 and 
^=7 t/ 4 (n=2) . (In (b) and (c) the angular variables were chosen so that the 
distortion of the wave was maximum.) As the figure illustrates, for a given 
value of M the amount of distortion decreases as n increases. Intuitively, 
we may think of this as being due to the increasing number of degrees of 
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freedom” of the wave. Alternatively, to obtain the same amount of distortion 
in the wave as n increases, the value of M would also have to be increased. 
This observation is consistent with our one-dimensional example (which can be 
thought of as a very "confining” flow), since for this example there was 
significant distortion of the wave at a much lower value of the Mach number. 

We also note that the analysis presented here may provide some insight 
into the derivation of appropriate non-reflecting boundary conditions to be 
used at artificial computational boundaries for a purely numerical simulation 
of acoustical waves generated by general (arbitrary) sources. For example, 
the rather simple form of our final results allows us to express the far field 
behavior of the waves in a concise form. This expression, when used with some 
of the ideas of Bayliss and Turkel (1980), for example, may allow us to derive 
the required boundary conditions. Investigations along these lines are 
continuing. 
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Appendix 


In this appendix we record the form of the perturbation coefficient 
functions for the velocity components, as well as various constants, 
associated with the examples presented in section 8. 

Pulsating sphere : 

u< 0) = (( r-1 )w/((l+o; 2 )r 2 )J cos (wr) + (l/(l+w 2 ) ) (w 2 /r + l/r 2 j sin(wr) , 

u^ = ( b_ / r 2 + 2b w/r)cos(2wr) + (b Q /r 2 - 2b w/r)sin(2wr) 

T Z o o jj 

+ ^(29-3'y)w 2 /(16r 3 (l+w 2 ) 2 ) + w 2 (w 2 -l)/(r 2 (l+w 2 ) 2 )j cos(2wr) 

- ( ( l+ 7 )w 2 log(r ) ) / ( 4r 2 ( 1+w 2 ) ) 

+ ^(29-37)w(w 2 -l)/(32r 3 (l+w 2 ) 2 ) - 2w 3 /(r 2 (l+w 2 ) 2 )j sin(2wr) 

+ 0(1/ r 4 ) . 

Oscillating sphere : 

uj 0) = ^2(r-l)w(2+rw 2 )cos(w7-)/r 3 + (4-2w 2 +4rw 2 -2r 2 w 2 +r 2 w 4 ) sin(wr ) /r 3 j 

• c o s ( 6 ) / ( 4+u 4 ) , 

u^°^= ^(w(2r-2-ru 2 )cos(wr)/r 3 + (2-w 2 +2rw 2 ) sin(wr) /r 3 j sin(0) /(4+w 4 ) , 

u^ 1 ^ R 0 (r)+R 1 (r)cos(2wr)+R 2 (r)sin(2wr) + 

^S 0 (r)+S 1 (r)cos(2wr)+S 2 (r)sin(2wr)j cos (20) , 

u^= ^T 0 (r)+T 1 (r)cos(2wr)+T 2 (r)sin(2wr)j sin(2^) , 
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f Qc = w 2 |2120+727+(U06-4307)w 2 +(519+77)w 4 +(408+247)w 6 ] /(384(1+4w 2 ) ) , 

f Qg = -u (^2028-207+( 344-6807)0 ; 2 - ( 1037-4997)u/ 4 

+(360+1047)w 6 -512w 8 j /(768(l+4w 2 ) ) , 

f 2c = w 2 ^-2106+1987-(141-5637)w 2 +(407-5537)o; 4 -(280+1047)w 6 

+ (408+247)w 8 j /(24(81+36w 2 -32w 4 +64o; 6 ) ) , 

f 2s = w(^4140-4687-(6136+2487)o; 2 -(4261-329l7)o; 4 + 

(96-4807)w 6 -(2512+3367)o; 8 +512w l0 j/(96(81+36w 2 -32w 4 +64u; 6 )) . 

Here each R ^ , Sj and T\ can be expressed as a multiple of (log(r))/r and a 
power series in the variable (1/r), beginning with a term proportional to 1/r. 

Squishing sphere : 

u£°^= sin(2f)sin 2 (0) |^( (4o; 5 -9u> 3 ) /r + (36o; 3 -4u; 5 ) /r 2 
+ (81w-36o; 3 ) /r 3 -(81o;-9w 3 ) /r 4 )cos(o;t) + 

( (81-36u; 2 ) / r 4 +(81u> 2 -9u> 4 ) / r 3 -(36w 2 -16u> 4 ) / r 2 -(9w 4 -u; 8 ) / r Isin^r) ) j 

/ (81+9w 2 -2o» 4 +o; 8 ) , 

u^^ = 2sin(2f)sin(0)cos(0) ^ 

( (27o>-3w 3 ) /r 4 -(27u;-12u; 3 ) /r 3 -(9u; 3 -u; 5 )/r 2 )cos(wr) + 

( (-27+12u/ 2 ) /r 4 +(9o; 2 -4u> 4 ) /r 2 -(27w 2 -3o> 4 ) /r 8 )sin(ow)j /(81+9o; 2 -2o; 4 +o; 6 ) , 

= 2cos(2^)sin(0) ^ 
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( ( 27w-3w 3 ) / r 4 - ( 27w-l 2w 3 ) / r 3 - ( 9w 3 -w 5 ) / r 2 ) c os ( wr ) + 

((-27+12w 2 )/r 4 +(9w 2 -4w 4 )/r 2 -(27w 2 -3w 4 )/r 3 )sin(u;r))/(81+9w 2 -2w 4 +w 6 ) , 


2 , 1 ^ 

p {2) = y~ , cos (2kg)|^cos(4j^) j+RC k .cos(2wr)+RS k>j sin(2wr)j j, 

i n ^ _ a 


k=0 \j=0 


= y cos(2kfl)j ^ cos(4j^) |s k j + ®^k j c 0 s ) +8^ j sin(2wr ) j j , 


k=0 j =0 


^ ‘ sin(2kfl)| ^ ' cos(4j^) j +TCj t j cos(2a;r)+TS k j sin(2wr) j|, 


k=l j=o 


z , 

u<, 1} = y~"sin((2k-l)g)|sin(4^)(u k 1 +UC fc 1 cos(2wr)+US k> 1 sin(2o;r) j 


k=l 


+ U k qH-UC, () cos(2wr)+US k 0 sin(2wr) 


Here each of the coefficients j, RC k j, ... , US k j can be expressed as a 

multiple of (log(r))/r and a power series in the variable (1/r), beginning 
with a term proportional to 1/r. 
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Figure 1 : Approximations to (p-l)/M for the simple one-dimensional 

example plotted as a function of distance from the source (at x=0) for 
t=42.44, using the multiple scales approximation from equation (4.18) (solid 
lines), the classical linear solution (4.15) (short dashed lines), and a 
solution to equations (4.1)-(4.2) obtained by purely numerical means 
(circles), with (a) M=0.01 and (b) M=0.04. 





Figure 2 : (a) An illustration of the Cartesian and spherical coordinate 

systems used in sections 5-8. Also, an indication of the vibrations of the 
sphere for the three examples considered in section 8: (b) the pulsating 
sphere; (c) the oscillating sphere; and (d) the "squishing" sphere, indicating 
the motion of the sphere in the plane 6- n/2. 


39 





























r (p~l ) /(Ms in 2 ( 0)s in(2^) ) t (p— 1 ) /(Msin 2 (0)sin(2^) ) 


Squishing Sphere 


-0.075 


$ = tt/12 


0.075 


^ = tt/6 


-0.075 


0.075 


$ = 7T/4 


-0.075 


0.075 


-0.075 


0.075 



w=l . 25 


t=50 

Squishing Sphere” 

M=< 


Figure 5 : Approximations to r(p-l)/(Msin 2 (0)sin(2^)) for the squishing 

sphere’' example plotted as a function of distance from the center of the 
sphere at U50 with w=1.25 and M=0.6, for (a) 0=xl4 and (b) 0=7r/2, at t = 

f 1 ) 2 

jr/12, tt/ 6, and *74, using the multiple scales solution rp /(3sin (0)sin(2^) 
(solid lines) and the classical linear solution (dashed lines). 







r(p-l)/M 


(a) Pulsating Sphere 


(b) Oscillating Sphere (6=0) 


(c) Squishing Sphere (0 = tt/2, V =ir/4) 


Figure 6 : Comparison of rp K k r(p-l)/M with M=0.3, t=50, and w=1.25 for 

(a) the pulsating sphere, (b) the oscillating sphere at 0=0; and (c) the 
squishing sphere at 6=i r/2 and ^=7r/4. In each case, the corresponding linear 
solution is plotted as a dashed line. 
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